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Abstract 

We study massive real scalar 0^ theory in the expanding Poincare patch of de Sitter space. 
We calculate the leading two loop infrared contribution to the two-point function in this theory. 
We do that for the massive fields both from the principal and complementary series. As can be 
expected at this order light fields from the complementary series show stronger infrared effects 
than the heavy fields from the principal one. For the principal series, unlike the complementary 
one, we can derive the kinetic equation from the system of Dyson-Schwinger equation, which 
allows us to sum up the leading infrared contributions from all loops. We find two peculiar 
solution of the kinetic equation. One of them describes stationary Gibbons-Hawking type 
distribution for the density per comoving volume. Another solution shows explosive (square 
root of the pole in finite proper time) growth of the particle number density per comoving 
volume. That signals the possibility of the destruction of the expanding Poincare patch even 
by the very massive fields. We conclude with the consideration of the infrared divergences in 
global de Sitter space and in its contracting Poincare patch. 



1 Introduction 



There are large infrared (IR) loop contributions even in the massive field theories in the 
Poincare patch (PP) of de Sitter (dS) space [I], [2], [3] (see also |1]-[S] for the massless fields). In 
global dS space there are IR divergences PP, [S], [ID] (see also [H]), which lead to the inevitable 
breaking of the dS isometry in the loops for any initial state. They are specific to dS space 
physics and are absent e.g. in AdS space [T2] . 

For the alternative point of view on the IR properties of the massive field theories in dS 
space see [13], [E]. However, these papers heavily relay on the analytical properties of the 
correlators as functions of the dS invariant distances between their arguments. But such an 
approach does not work when the dS isometry is broken by the initial density perturbation (or 
by the IR divergence in global dS). In fact, then the correlation functions stop to depend only 
on the invariant distance and become functions of each point separately. The motivation to 
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include initial density perturbation is that it is rather unphysical to address the question of the 
stability of the system under such conditions when all its symmetries are respected exactly. We 
propose to consider excitations above the highly symmetric state and to trace where do they 
evolve in the future infinity — to the symmetric state or not? In global dS one even does not 
have to consider initial density perturbations, because dS isometry is broken in the loops by 
the IR divergences. 

In m, [2], [3], [9], [in], [IS] the large IR contributions to the two point functions were 
calculated for the massive real scalar field theory with the cubic, 0^, selfinteraction. The cubic 
potential has the runaway instability. Moreover, one-loop higher power correlation functions 
of the 0^ theory may potentially show strong IR effects, which are absent in the 0^ theory 
at least in 4D. To show that those IR effects, which are observed in our previous papers, are 
universal rather than due to do the specific properties of the 0^ theory, we consider here the 
scalar field theory with 0^ selfinteraction. Furthermore, in the pervious papers only massive 
fields from the principal series have been considered. Here we extend those considerations to 
the complementary series. 

We calculate loop corrections to the so called Keldysh propagator in PP of dS space. Unlike 
0^ case, in the 0^ theory there are no any large IR contributions in the first loop. However, at 
the two loop order we find such contributions in the sunset diagrams. In the case of the fields 
from the principal series the contribution is linear logarithmic in the physical momentum. For 
the complementary series the situation is more subtle and depends on the type of a-harmonics, 
but always is power like, i.e. stronger than for the principal series. 

We are not yet capable to perform the summation of the higher loop contributions for the 
case of the complementary series. But for the principal series it is possible to do such a summa- 
tion. That is done via a suitable IR ansats for the solution of the system of Dyson- Schwinger 
(DS) equations, which allows to reduce it to a generalization of the Boltzmann's kinetic equa- 
tion. The latter one has a clear physical meaning and describes various particle decay and 
creation processes in the dS space (see also [I7|-[20], [H]). We solve this kinetic equation 
in two cases. One corresponds to the very mild initial density perturbation over the initial 
Bunch-Davies state. In such a case the state of the theory relaxes to the eventual Gibbons- 
Hawking type stationary distribution for the number density per comoving volume. Another 
situation corresponds to the strong enough initial density perturbation, which, however, is still 
much smaller than the cosmo logical constant. In this case the state of the theory shows the 
explosive (square root of the pole in finite proper time) growth of the particle number density 
per comoving volume. 

We conclude with the consideration of the contracting PP of dS space and of the global dS 
and draw similar conclusions to those which have been made in [lOj for the 0^ theory. 

2 The setup of the problem 

L^-dimensional de Sitter (dS) space-time is the hyperboloid, Xq — Xf = —1, i = 1, . . . ,D, 
inside [D + l)-dimensional Minkowski space-time, ds"^ = dX^ — dXf. Throughout this paper 
we set its curvature to 1 and mostly consider its half (e.g. Xq > X^), which is referred to 
as expanding Poincare patch (PP): ds^ = ^{d^f — dx"^), where rj G (-|-cxd,0). Note that while 
7/ — > +CXD is the past, 77 — )■ corresponds to the future infinity. The action of the theory which 
we are going to study is 
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S= I d^x 



Throughout this paper we always assume that m > 0. 

Corresponding free harmonics in the Fourier expansion 77) = J d^~^krj^~ \ak h{kvi) e~^^^ 

+al h*{kri) e*^^] , k = k are defined via h{x), which is a solution of the Bessel equation with 



the index i jj, = i \J m? — (-^^) • The choice of such a solution specifies the dS invariant 
vacuum \vac) = 0. E.g. Bunch-Davies (BD) state (or in-vacuum in PP) corresponds to 
h{pri) = ^ H^^\pri), where H^^^ is the Hankel function of the first kind. All other 
harmonics and vacua can be obtained from those of BD via a one parameter family of the 
Bogolyubov rotations [21], [22]. For example, the out- vacuum corresponds to the so called out- 



Jost-harmonics, h{pri) = y -^^^^ Jifi{pv) j where Ji^ is the Bessel function of the first kind. 

Conjugate harmonics in the latter case are given by the F's — Bessel functions of the second 
kind. 

In the nonstationary situation every particle is described by the matrix propagator (see 
e.g. [23] for the Feynman rules in the 0^ scalar field theory on dS space) whose entries are — 
the Keldysh propagator, Gk = \ ({0(3^), 0(y)}), and the retarded and advanced propagator^], 

= =F([0(a;), (p{y)]) 6'(=F^a^o)- The more detailed discussion of the physics in dS space, which 
is relevant for the present paper, can be found in [1], [3], [9], [10], [15]. It is instructive for 
the further discussion to keep in mind that the result of the calculation with the use of the 
non- stationary Keldysh-Schwinger diagrammatic technique provides a solution to a sort of the 
Cauchy problem. 

Because PP is spatially homogeneous and because below we concentrate on the spatially 
homogeneous states, we find it convenient to make the spatial Fourier transform of the afore- 
mentioned propagators: 

D''''''^ {Vi,V2,p) = I d^-'xe'P^'G^''''^ {Vux,V2,0), 

D^{Vi,V2,p) = {ViV2)^d^{prii,pr]2), 
Di{vuV2,p) = T6{±^v){ViV2)~d-{prii,pri2), (2) 

where Ar] = rji — ri2, p = \p\ and 



d {pr]i,PV2) = 2Im [h{pr]i)h*{pr]2)] , 
d^{pr]i,pr]2) = h{pr]i)h*{pr]2) Q + rij^ + h{pr]i)h{pr]2)kp + c.c. . (3) 

Here we use the convenient notations = (a^ Op) and Hp = {ap a_p) . In the non-stationary 
situation Up and kp can be zero only on the tree-level, i.e. if the initial state is chosen to be 
\vac). 

It is worth stressing here that in our paper we always study two-point correlation function. 
According to ([2])-([3]), the quantities Up = {ap ap) and kp = (apa^p) are juts elements of this 

*Here Aa;o = — j/Oi {7 } is the anti-commutator and [, ] is the commutator; 9{ri) is the Heaviside ^-function. 
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correlation function. Hence, all our conclusions about their behavior have invariant physical 
meaning. Furthermore, in the situation when the anomalous quantum average kp vanishes Up 
acquires the clear physical meaning — it becomes the particle number density per comoving 
volume in the given state of the theory. In fact, then free Hamiltonian of the theory acquires 
the diagonal form. Moreover, then all elements of the kinetic equation which is presented below 
acquire the clear physical meaning. 

Let us provide here a few more arguments favoring our interpretation of np as the particle 
density. Various observers may indeed detect different particle fluxes. However, one should sep- 
arate Unruh effect from what we would like to call as the real particle production. In Minkowski 
space both inertial and non-inertial observers see the same state - Minkowski (Poincare invari- 
ant) vacuum. However, while inertial observer sees it as the empty space, the non-inertial one 
sees it as the thermal state. That is due to the specific correlation of the vacuum fluctua- 
tions along its world-line [I^. Note that there is no any non-trivial gravitational field in the 
circumstances under consideration, because Riemanian tensor is exactly zero. 

The real particle creation is due to the change of the "vacuum" state under the influence 
of the background field. That is exactly what happens in the background electric field, in dS 
space and in the collapsing black hole background. 

Rephrasing, we would like to say here that while in Minkowski space there is one type of 
observers which does not see any particle flux, in dS space there is no such an observer which 
sees nothing. On general grounds we expect that the least particle flux is seen by inertial 
observers — they do not see extra Unruh type of flux, so to say. Our calculations are actually 
done for the inertial observers, because the change from the proper time t to the conformal 
one, 7] = e~*, is just the change of the clock's rate rather than a transition to some non-trivial 
motion. In any case all our formulas can be trivially rewritten in the proper time t. 

Finally, the IR dynamics depends on the value of the mass m and on the choice of h{x). 
Below we separately consider the following two cases: ('^^)^ < — the principal series, and 

{^j^Y ^ — complementary series. The crucial physical difference between these two 
cases is due to the fact that while the harmonics of the principal series oscillate at the future 
infinity, h{x) ~ A+x*^ + A_ x^*'^, x 0, those of the complementary series don't do that, 
because in the latter case fi is pure imaginary. 

3 Two-loop contribution 

In this section we calculate loop corrections to the Keldysh propagator, {r]i,r]2,p), with 
the initial dS invariant vacuum state at the past infinity of PP. In these settings the tree-level 

{rii,ri2,p) is given by ([3]) with Up = and Hp = 0. It is straightforward to show that 
the 0^ theory, unlike 0^ one, does not poses any large IR contributions to any propagator 
at the first loop order (~ A). However, in the second loop order (~ A^) there is a large IR 
contribution to , which is of interest for us. We consider the IR limit in which p^/rjirfi 
and rii/rj2 = const [3], [15], [10]. It corresponds to the situation when both time arguments of 
the propagator are taken to the future infinity, while the time distance between them is kept 
finite. 

The reason why we pay attention only to the Keldysh propagator is that it defines the state 
of the theory, i.e. shows the dependence of Hp and kp on time. Moreover, in the case of principal 
series it is the only propagator which receives large corrections of the order A^ log(p ^771772/ /^) in 
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the IR limit in question. In fact, in [2J it was shown that the retarded and advanced propagators 
of the 0^ theory receive only finite corrections from the first loop — of the order log(?7i/?72). 
In the case of the 0^ theory the situation is similar in the second loop. Furthermore it is 
straightforward to show that for the principal series the large IR contributions to the interaction 
vertex are also suppressed by higher powers of A. For the complementary series the situation 
is more subtle, but we still would like to start their study with the consideration of the loop 
corrections to . 

The large IR loop contribution to (if there is one) at the order A^ can be expressed in 
the form ([2]), ^ with 



A^ f n n 

^piv) - ~ 3(27r)^(^-i) J c?^~^g3 J d7]3 J dr]i{7]^r]i)^-'^ 

5^^'^^ (p + gl + g2 + gs) h{pri:i)h{qiri^)h{q2ri:i)h{q^ri^)h* {prii)h* {qivji^^^ 

kpiv) ^ 3 (2^)2(D-i) J d^^^'q.d^-^d^-^ j dvsj dr^.ir^sV^)^-' 

^(^5-1) (^p- _^ _^ ^2 + (fs) h%pr]3)h{qir]3)h{q2r]3)h{qsr]s)h%pr]4)h%qir]4)h%q2r]4)h*{q3r]4). (4) 

Here j? = gi,2,3 = |'fi,2,3| and np{ri) is real. Note that if one will take the flat space limit of 
these expressions, i.e. substitute t] by t, \/\g\ by 1 and r/('^~i)/2 h{k7]) by e"*^'^)*, then Up and kp 
would vanish, when — )■ 0, as the consequence of the energy conservation. 

In deriving these expressions we have used that h{pr]) depends only on |p| and, hence, we 
can safely change p —p. Also in the IR limit in question inside the leading loop corrections 
one can neglect the difference between rji, 772 and 77 = y/rji 772- In such an approximation we drop 
the subleading, A^ log {rii/ri2), contributions from the expressions for Up and kp. The derivation 
of (jl]) is similar to the one performed in P, [2], [3], ^0\, [15] for the (j)^ theory. 

3.1 Principal series 

We start with the case < m and make the following change of the integration variables 
in (jl]): (ji to li = (jiTj^ and 774 to f = ^, i = 1, 2, 3. Then we expand hijyq^^/i) ^ ^+(^^73,4)*'^ + 
A_{pri3^4)~^^ as pr73,4 — )■ under the integrals, where A± are some mass dependent complex 
constants. After that we neg lectll] p in comparison with gi 2,3 on the RHS of (jl]) and perform 
the integration over rj^. The largest IR contributions come from the integrals of the type 
//Vp {PV3) h* {p^zv) and JJ^^ (PVs) ^* {PVs'^)^ where /I's are Taylor expanded. The 
result is as follows: 



Up, ^ / dvv^-' I d^'-'h d^'% d^-% (/-; + + h) X 

xh*{h)h ihv) h*{l2)h ikv) h*{h)h ikv) [\A+f v-'^ + \A^f v^'^'] , 

kp, ^ ^^(2^^!)^ /' dvv^^' I d^-% d^-H2 d^-% (/-; + + h) X 



xh*ih)hihv)h*il2)hil2v)h*{h)hihv)A+A_[v'^ + V-''^], /i,2,3= ri,2,3 • (5) 



^The justification of all the listed here approximations can be found in [T], [2], [3], [5], [TS], [lOj . 
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The lower limit of integration over 773 we cut by yU, because at prj ^ jj, the d drj^ and 
dv integrals are rapidly convergent due to the oscillations of h{x), while we care only about 
the leading IR contribution. It is worth stressing here that such a convergence of the dv and 
d^~^li integrals is true in the sense of the generalized functions. The latter fact is related to 
the discussed below behavior of the modes with high momenta. 

Note that 77,^(77) and kp{'q) are functions of the physical momentum, pr], rather than sepa- 
rately depend on the momentum, p, and on the time, 77. That is natural because of the spatial 
homogeneity and the invariance of the PP dS metric under the simultaneous rescaling x — )■ ax, 
r] — )■ at]. 

It was explained in |3], [9], [10], [15] that such large log contributions in appear due 
to the particle creation in dS space (see also the discussion below). In the expanding PP the 
creation of particles with comoving momentum p effectively starts after some moment of time 
^* ~ /i/p, because the modes with high momenta, prj ^ /i, do not feel the curvature of the 
space-time and behave as if they are in fiat space. As the result Up and kp are proportional 
to the proper time, log(/i/p?7), elapsed from [x/p to rj. The coefficient of the proportionality 
is just the particle production rate. It is worth stressing here that the presence of the large 
kp, which is comparable to rip, signals that there is the strong backreaction on the initial state 
\vac) ^MM- 

For the further reference it is instructive to study also the IR behavior of Up and kp for 
the out-Jost-harmonics. Although UV behavior of the correlation functions for the out-Jost- 
harmonics is different from the proper fiat space type, we are interested here in the IR limit, 
where out-Jost-harmonics may play a crucial rolqil. In fact, in condensed matter physics it is 
the frequent situation that in the IR limit one has to perform a Bogolyubov rotation to some 
harmonics whose UV properties may be different from the proper ones, however they properly 
describe the IR physics. The seminal example is the BCS theory for superconductivity. 

In particular, we will see that the proper solution of the IR limit of the system of Dyson- 
Schwinger (DS) equations is obtained via the out-Jost-harmonics. The hint for that comes from 
the following observation. For the out-Jost-harmonics the leading IR two-loop contribution has 
a crucial difference wrt that of BD harmonics or any other a-harmonics. In fact, the out-Jost- 
harmonics behave as h{p7]s^4) ^ ^(^^73,4)*^ in the future infinity, prjs^A — )■ 0, where A is some 
mass dependent complex constant. Then it is straightforward to show that Up has the same 
form as ([5]) with instead of [|v4+|^ t>*'^ + \A^f v~^^'\. At the same time kp does not 

receive any large contributions in the IR limit in question, i.e. it is negligible in comparison 
with Up. This is going to be an important observation for the derivation of the Icinetic equation 
below. 

It is worth stressing here that for the out-Jost-harmonics the dv and d^~^li integrals in the 
obtained expressions for Up and kp are also convergent. In fact, the situation in the 0^ theory 
is even simpler than in the 0^ case [3]. 

3.2 Complementary series 

We continue with the consideration of the complementary series, > m, corresponding 
to the imaginary /z, i.e. to the real index of the solution of the Bessel equation h(jyq). Below 

■fWhile the theory under consideration shows the proper UV behavior in the BD state, it does not do that 
in any other a-state. The reason for that is as follows. While BD harmonics behave as single waves e'^'' in the 
UV limit, prj — > 00, the other a-harmonics behave as linear combinations of e'^'' and e"'^'' in the same limit. 
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we use the notation v = —i^. Then for the in-harmonics h = + lYy, where both Bessel 
functions Jp and Yjj are real in the case of the real index v. Expanding them near zero, we get 
Y'^(x) ~ A^x^'^ + Bx^^^"^ and Jy{x) ~ A^x^ . Because of the possible differences between the 
behavior of h and h* near zero we have to pay attention separately to kp and fc*. 
The contributions to and fcp, k* can be expressed as 



kp(r]) 



2A2 



du 



3(2 7r)2(^-i) u 



du 



dv 



.2D' 



u 



jF[v\h{uv)h* [-) e[uv-pr]]e 



dv 



,2D 



Kiv) 



4A2 r^du r° dv 



— F[v]h*{uv)h* 1-] e[uv-p7]]e 



u 



3(27r)2(^-i) u 



,2D- 



j F[v] h{uv)h i — j 6[uv — prj] 6 



u 



prj 



u 

uv 

V 

~u 

uv 

IV 



(6) 



with the use of the following notations 



M4)^-' 6^"''^ {p- qi - q2 - gs) , (7) 

which, after the change of integration variables u = p^/rj^rjl, v = can be expressed as 

F{r)3,r)4) = ^^^-^2(0^1) M^^'^"^-' F[v], where F[v] = F* [l/v] is some function of one variable v. 

The leading correction to Up and kp is given by ([6]) where from the Hankel functions, h{uv) 
and h{u/v), we single out only F's. Such a contribution gives for Up and kp the inverse power 
like behavior in prj, which, however, cancels out after the substitution into , because Y is 
real. The next order is obtained as follows. One also has to express /^(pr/i 2) through J^, and 
Yi, in the full propagator . Then from one of the four /I's {h{prji^2) and h{uv), h{u/v)) we 
have to single out J^,, while from the other three — Y^^s. This expression doesn't cancel out 
and provides the leading IR contribution to D^: 



Y\_d^ ^QMqi'n3)h*{qir]i) 



i=l 



D {vuV2,p) 



3(27r)2(^-i) 



du u 



-l-2u 



dvv^-'^'Flv 



-1 + 









/ . A 


e 


w 




uv 


yprjv J 




V. 



9 [—prj + uv] + 



2v 



U 

-UV H — 



u 

-prj + - 

V 



After the straightforward manipulations the obtained expression can be reduced to: 



D^iVi,V2,p) 



8 A^ A^_ A+ rj^~^ log(p?7) 
3(27r)2(^-i) (pry)2- 



dvv 



1-2D 



F[v] 



2u \v 



2v 



1 

- jdw'-''' 





F[v] 



— (vY 
2u ^ ' 



2u ^ ^2. 



(9) 
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The integral over v is convergent in the IR hmit (as v — > oo) if D > 1 + 4i/. In the UV hmit 
{v — > 0) it is convergent in the sense of the generahzed function. 

For the out-harmonics the situation is a bit different. In this case h — J^,, h* — Y^,. The 
straightforward calculation shows that 



■2V+1-2D 



np{ri) oc log{pri) / dv F[ 

Jo 

kp(r]) oc A' At (pr))-^" J dv F[v] v^ 

k;{r))o^X^Alii^'' [' dvF[v]v-^''+^-^'' + 0{{pr)f''). (10) 
Jo 

After the substitution into the Keldysh propagator the leading contribution comes from k* and 
is as follows: 



\2 A2 A2 2u pi 

i;i^~_^J2^_li±^i^^^-i / dvv-^^+^-^'^Fiv] (11) 
3i/(27r)^-i(pr;)2- ' ^ ^ ^ ^ 



Thus, for the complementary series the loop contributions are power like, i.e. stronger than for 
the case of the principal series. 

Because of the character of these IR contributions for the light fields we do not yet under- 
stand their physical meaning and think that the obtained in the next section kinetic equation is 
not applicable for the fields from the complementary series. In the case of the complementary 
series we do not yet know how to perform the summation of the leading IR contributions from 
all loops. But on general physical grounds and from the two loop result we expect that com- 
plementary series will show stronger IR effects than the heavy fields from the principal series, 
which is under study below. 



4 Kinetic equation 

Although A^ is small the product A^log {fJ^/pi]) can become large as prj — > 0. Hence, higher 
loops are not suppressed in comparison with the tree-level contribution. Then, one has to 
perform the summation of the leading IR contributions from all loops. This can be achieved via 
an approximate IR solution of the system of DS equations in the non-stationary diagrammatic 
technique. We can do that only for the case of the principal series, m > (D — l)/2, at least 
because in this case harmonics oscillate in the IR limit. As the result there is a clear separation 
of scales between the time dependence of the harmonics h{pri) and that of and kp^. This 
allows to simplify DS equations in the limit under consideration and eventually to solve them 
approximately. In the case of the complementary series, however, the character of the IR 
contributions to and to the vertex does not allow to have a clear kinetic interpretation. 

Putting it in other words, to move further it is worth observing that above we have calculated 
the loop corrections to the Keldysh propagator under the assumption that Up and kp retain 
their initial values throughout all the time evolution. In fact, in the calculations of the previous 
section we have used tree-level propagators, i.e. their initial values. To make the problem 
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self-consistent one has to take into account the change of and kp in time. As we will 
see, this also allows to reduce the problem to the solution of the DS equations for the non- 
stationary diagrammatic technique. The latter ones represent a system of equations for the 
matrix propagators, self-energy and for the vertex. In some circumstances, as we will see, this 
system can be simplified and reduced to a single equation. 

We assume that the evolution had started with some density perturbation over the BD 
state at the past infinity of PP. Then the UV behavior of the theory in question is the same 
as in flat space, but dS invariance is slightly broken. We would like to trace the destiny of 
these density perturbations in the future infinity, i.e. would like to see whether the theory 
relaxes back to the dS invariant state or this density explodes causing the modification of the 
background geometry? The answer to the latter question also can be obtained only after the 
solution of the IR limit of the non-stationary DS equation. 

Thus, we would like to sum up leading contributions, which are powers of A ^ log {py/f]i V2/ f^), 
and drop the sub leading terms, which are suppressed by higher powers of A and/or powers of 
A^ log (771/772)- Having in mind that retarded and advanced propagators and the vertex do not 
receive any large contributions in the limit in questioiiH, we can assume that they take their 
classical (UV renormalized) values. At the same time we assume that the ansatz for the exact 
Keldysh propagator is given by ([2]), ([3D with undefined np{r]) and kp{ri), where r] = ^7/1772. 

Then the relevant part of the system of DS equations has the form 



A2 



D-l. 



D^{r]i,r]2,p) = DQ{r]i,r]2,p) - 
-5^ >{p-qi-q2- ga) 



3 D^iVu V3,P)D^{V3, V4, qi)D^{r]3, 7/4, q2)D^{r]3, 774, q3)D^{T]^, 7/2, p) 
1 



-^DQ{r]i, r/3, p)Do{rj3, 774, gi)L)^(r/3, 774, g2)^^(??3, Vi, ?3)^o^(^4, V2,p) 
3 

D§iVuV3^P)D^iV3, V4, qi)D^{7]s, r/4, q2)D^{r]3, r/4, g3)D^(r/4, r]2,p) 

+Do{vu V3, P)D^{V3, V4, qi)D^{r]3, 774, 52)^^(^/3, V4, q?i)D^{Tli, V2,p) 
3 

DoiVl, V3,P)D^iV3, V4, V4, q2)D^iV3, Vi, 93)^^(^4, V2,P) 



-^^0 iVi^V3, p)Dq (r/3,774,gi)£'o (^3,^4,52)^0 iV3,Vi,(l3)D ipVi,PV2) 



+3 Dl^{r]i, r/3, p)D^{r]3, r/4, gi)D^ (773, 774, ^2)^"" (^3, ^4, 93)^'' (^4, V2,p) 



K I 



(12) 



where p = \p\, ^1,2,3 = |'?i,2,3| and Dq'^'^ are the Fourier transforms of the initial values of 
the retarded, advanced and Keldysh propagator^, — is the exact Keldysh propagator. 

''Of course there are large IR contributions to D"^'^ and to the vertex, which are coming from those in 
. Correspondingly they are suppressed by the higher powers of A, because they appear in the higher 
loops. Actually all one-loop contributions to the vertex have similar structure to the one shown in with one 
momentum integration less. It is straightforward to see that because of that there is no any large IR contribution 
to the vertex at the order. 



^I.e. is also given by ([I]),® with some initial values n^"^ and k 



AO) 
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This equation is covariant under the Bogolyubov rotations between different harmonics, h{x). 
Because we are interested in its solution in the IR hmit we do not have to care about the proper 
UV behavior and have to check the situation for all possible a-harmonics. 

We would like to pick out the largest IR contribution from the integral on the RHS of f|T2|) . 
The calculation is just a straightforward generalization of the above two-loop one. Similarly to 
j3], [10], |15], for all «-harmonics, including BD ones, we obtain that the ansatz (j2]), (j3]) solves 
DS equation in the IR limit in question. But kp is comparable to Up. This means that the 
backreaction on the background state (specified by the choice of the harmonics) is big. The 
only exception is given by the out-Jost-harmonics. For these harmonics kp remains zero if its 
initial value was zero. Moreover, small perturbations of kp relax back to zero. For the 0^ theory 
this was shown in [15] . For the 0^ theory the situation is similar. 

As the result, with the use of the out-Jost-harmonics the IR limit of the DS equation is 
solved by (EJ), ([3]) with kp = 0. We would like to convert the integral DS equation (JT2ll into the 
integrodifferential form, i.e. into the form of the kinetic equation [21]. This is done via a kind 
of the renormalization group procedure as follows [15] . In the given settings n^^ is the particle 
density at some moment after r^^, ~ n/p. In fact, as we have mentioned above and explain in 
the next section, before this moment all the kinetic processes, except may be the irrelevant 
scattering one, are suppressed. Hence, np{ri) remains constant and is equal to Then, 
from f|T2|) it is straightforward to derive that the difference between np{ri) and n^p^ = np{r]^) is 
proportional to the proper time elapsed from r/* to rj. The coefficient of the proportionality is 
the collision integral — the RHS of the kinetic equation: 



dn. 



PV 



np{r]) - np{r]^) 

log(?7) - log(?7=,) ' d\og{pri) 

v'" h*{k)h*{l2)h 



3Re 

X 

-|-3Re 

X 

+Re 

X 

+Re 



X^\A\' 
6 

h + k 



D-l 



(27r)^-i (27r; 
h{hv)h{l2v)h* 



k + k 



(1 + np^)ni^ni^{l + n^r^^f^^) - ^^^(1 + ri/J(l + ni^)n^^^^^^^ 



v'^ h*{k)hik)h 



k-k 



h{kv)h*{kv)h* 



k-k 



(1 + np^)ni^{l + n;J(l + n^i^_^^^) - n^^il + ni^)ni^n^^^_^^^ 



k + k 



h{kv)h{kv)h 



k + k 





.,)] 


X 


^\h+l2\_ 




X 




r2i_ 




X 



v'^ h{k)h{k)h 



k + k 



h*{kv)h*{kv)h* 



k + k 



X 



X 



(13) 



In the process of the derivation of this equation we have neglected p in comparison with qi^2,3 
on its RHS, denoted k,2,3 = (li,2,3V assumed that on the RHS of f|T2|) rifc^ is much slower 
function of time than h{kri). The latter fact is true because of the above mentioned separation 
of scales. As the result we can safely take out all n's from the argument of the time integral on 
the RHS of ([H]) and substitute by . 

Note that one can reproduce Hp in eq. ([5]) from ( !T3|) if he will use Hankel functions in place 
of h{x), change lAjV' to [\A+\'^v'f + lA^l^V'^"] and put all n's to zero on the RHS of ([13]). 
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5 Solution of the kinetic equation 



The kinetic equation (fT3|) does not posess Plankian distribution as its solution, because of 
the violation of energy conservation in the dS background gravitational field. It can be mapped 
to the kinetic equation in flat space via the substitution of the harmonics t^^-^"^)/^ h{kri) by 
the plane waves e*'^'-'^-** in the collision integral. In the latter case one has 5-functions ensuring 
energy conservation on the RHS of ( IT3l) instead of the integrals of /I's. That in particular means 
that for high energy modes, krj ^ /i, the kinetics in dS space is the same as in fiat one: they 
can scatter off each other, but cannot be created in various processes involving dS background. 

Thus, suppose we have started at past infinity of PP with some very mild density perturba- 
tion over the BD state. After the Bogolyubov rotation to the out-Jost-harmonics, one has some 
initial values of np and kp. As can be understood from the discussion in the previous paragraph, 
for the given p the density Up and the anomalous average kp practically do not change before 
the moment r^* ~ fi/p. After this moment they begin to evolve according to the coupled system 
of kinetic equations for Up and kp. (To simplify the presentation we do not show them here. 
Similar equations for the (p^ theory can be found in [3], [IQ].) If kp is sufficiently small it relaxes 
to zero [T5] and the problem is reduced to the solution of (fT3ll . 

Now if the initial value of Hp after the rotation to the out-Jost-harmonics is much smaller 
than 1 we can use the following approximations: 



(1 + np^)ni^ni^{l + n^r,+h\) " ^mi^ + + ni^)n'\r,+r2\ ^ 

(1 + npr,)ni,{l + n/J(l + r^|fJ_r2|) " ^pvC^ + ^11)^12^1^-12] ~ 
(1 + np^)ni,ni^n^j;^_^_i'^^ - np^{l + + + n^f^^i'^^) ~ -n^^ 

(1 + nprj){l + + n,J(l + n^ri+i2\) " ^Pv^hnhnin+r^l - 1- 



(14) 



Because of the rapid oscillations of h{x) as x — 00 the integrals on the RHS of ( IT3|) are 
saturated at li ~ fi. At the same time prj <^ fi. Furthermore it is natural to assume that 
n/^^i <^ Upri in the situation of the very small initial density perturbation and the vanishing 
production of the high momentum modes. Hence, we can neglect the second term on the RHS 
of (1131) in comparison with the third and the fourth ones. 
As the result (ITS]) is reduced to: 



dn. 



d\og{pri) 



6 



(27r)^-i (27r) 



D-l 



xRe 



xRe 



v"" h*ih)h*{l2)h* 



h{kv)h{kv)h 
2 



D-l 



v'" h{k)h{k)h 



k + k 



(27r)^-i (27r; 
h*{kv)h*{kv)h* 



T2 , where 







dv V 

00 

k + k 


dv V 
k + k 



D~2 



X 



D-2 



X 



(15) 



Here Fi and F2 are the particle decay and production rates, correspondingly. Note that prj is 
reducing to zero in the approach towards the future infinity. 
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The obtained equation (fTSl) has the solution with the flat stationary point distribution 
npr, = r2/ri, which corresponds to the situation when the production (gain) of particles on 
the level pr] is equilibrated by the particle decay (loss) from the same level. Note that here we 
are talking about the number density per comoving volume, which would stay constant if there 
were no particle decay and production processes. 

The obtained solution is self-consistent for the large enough /i, because then r2/ri ^ 
g-3 7r/i <^ 1. Note that the equilibrium distribution is n g-STr^t ^^^^ jg quite Gibbons- 
Hawking one. Apparently, it looks like the thermal Boltzmann one, but the temperature 
depends on the power of the selfinteraction potential. In fact, the stationary distribution in the 
0^ theory is n ~ e"^'^'^. 

So the result of the summation of the large IR contributions may lead to the finite exact 
two point functions. But what if the evolution had started with some quite strong density 
perturbation (which is, however, still smaller than the cosmological constant) over the BD 
vacuum state? Now we are going to show that there is another very peculiar solution of the 
kinetic equation under consideration. See [3], [10] for the similar discussion in the case of (p^ 
theory. 

Suppose that due to the particle creation by the background gravitational field and by the 
particle decays from the other levels the density per comoving volume on the given level with 
pr/ -C /X became big in comparison with 1. Taking into account the flatness of the spectrum 
in dS space, it is natural to expect that for the harmonics with the low physical momenta the 
density very slowly depends on its argument. Hence, we can assume that n{pr]) ^ n{qi^2,3V) 
prj <^ ^ and gi, 2,3^7 <^ /U. 

Then, we can make the following approximations: 

(1 + npn)ni^ni^{l + ni.^) - ?2p^(l + n/J(l + n^Jn^g ^ 
(1 + rap^)n,,(l + n/J(l + ni^) - np^{l + ni^)ni^ni^ ^ 2ra^^ 
{I + np,^)ni^ni^ni.^ -np^(l + n;J(l + ni^){l + ni^) ^ -2nJ^ 
(1 + np^)(l + ni,){l + ni^){l + m^) - npr^ni^m^ni^ ^ Anl^ (16) 

and accept that on the RHS of ( !T3|) the main contribution to the U integrals comes form the 
region where U <^ /x, because n{x) ^ n{ii) ii x <^ n and y ^ ji. As the result the kinetic 
equation reduces to: 
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dn. 



vv 



d \og{pri) 



-TUp^, where 



3 UIMRe 



(27r)^-i 

ifi 



(27r)^-i 



(if f 



_D-2 



h{hv)h*{l2v)h* 



iA*f Re 



+2 A-^ Re 



k + h 



h+k 



h{kv)h{kv)h 
h*{kv)h*{kv)h 



k-k 

k + k 
k + k 



Note that T is independent of p. This equation has the solution: 



(17) 



ripv ~ , - (18) 

Q 

where 77* = ^e~2r and C is the integration constant, which depends on the initial conditions. 
The obtained solution is valid if /x/p = > > r/^. 

Thus, we see that there is a singular solution of the kinetic equation under consideration, 
which corresponds to the explosion of the particle number density per comoving volume within 
a finite proper time. Of course such an explosion wins against the expansion of the PP, because 
jjK ^ fj^-^ I y/\ogri/ri^,. Hence, the energy-momentum tensor of the created particles becomes 
huge and the backreaction has to be taken into account. As the result the dS space gets 
modified. But that is the problem for a separate study. At this point we just would like to 
stress that we see catastrophic IR effects even for the massive fields. It is natural to expect 
that for the light fields the situation will be even more drammatic. 

6 Comments on the contracting PP and global dS (in- 
stead of Conclusions) 

Contracting PP of dS space is interesting at least because it is complementary to the ex- 
panding PP within global dS space and we find it quite dangerous to study such a geodesi- 
caly incomplete subspace as PP alone. Contracting PP is represented by the same metric, 
ds"^ = ^ [dr]"^ — dx^], as the expanding one, but now the conformal time r] flows in the proper 
direction — from zero, at the past, to infinity, in the future. 

From f|T3|) one can straightforwardly obtain the kinetic equation in the contracting PP, 
if he will consider perfectly spatially homogeneous states. The latter situation while being 
stable in the expanding PP, is unstable in the contracting one under small inhomogeneous 
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density perturbations. However, it is still instructive to consider such an ideal situation in the 
contractog Pli, In this section we restrict ourselves to the case of the principal series. 

To perform the map between the expanding and contracting PPs both in the few loop 
calculations and in the kinetic equation one just has to flip the limits of drj integrations and 
change h{x) and h*{x), because of the exchange between positive and negative energy states 
under the flip of time. 

Then, it is straightforward to see that loop corrections to have the explicit IR divergence 
because particle creation starts right at the moment r^o, when we switch on selfinteractions. The 
divergence reveals itself via the impossibility to move the r]o to the past infinity. Due to the 
blue shifting of all modes, however, all the relevant kinetic processes stop after r]^: ~ fi/p. As 
the result the two loop divergence in question is proportional to log(?7/?7o)5 when prj fi and to 
log{fi / prjo) , when prj ^ fi. The prefactors are easily derivable, in view of the above discussion. 
Also it is straightforward to show that for the in- Jost- harmonics (Bessel functions in place of 
/i's) of the contracting PP the kp behaves similarly to that of the out-Jost-harmonics of the 
expanding PP. 

In conclusion, in the contracting PP we can find similar solutions of the kinetic equation to 
those, which have been found in the previous section. E.g. eq. (fT^ is mapped to 



-^^^^-r,n(r^) + r2, (19) 
dlog{r]/r]o) 

which shows a very peculiar phenomenon that despite that in two loops rjo — the moment of 
switching on selfinteractions — cannot be taken to the past infinity, after the summation of all 
loops we may find the theory at the stationary point state, Up = r2/ri, which allows one to 
remove the IR cutoff, r]o. 

At the same time the solution (fT8l) is mapped to: 



np{7]) ^ _ =, (20) 

V2r log {r]^/r]) 



where 77^ = 7706"-^/^^ < yu/p. Of course whether the field theory state goes into ( |20l) or to (fT9l) 
depends on the initial conditions. 

The situation in the global dS space is even more interesting. On the one hand all sorts of 
contributions to the collision integral in global dS space are sums of those in expanding and 
contracting Pp£]- But depending on the choice of harmonics we either have IR divergence 
in kp, or kp is divergent as the system advances towards the future infinity. Or even more 
generic situation for a- vacua is when kp has both types of such divergences simultaneously. As 
the result there is no any choice of harmonics in which kp is negligible in comparison with Hp, 
which probably means that there is no any stationary state in global dS and the backreaction 
on the background geometry is strong with any initial conditions. 



"Actually it is not very hard to find the inhomogeneous extension of the kinetic equation ([T^ . In the case 
when the particle density starts to depend on the spatial position Up — np{x) one has to substitute rjd/dri on 
the LHS of (dH) by ?7 9^ + r]^pd^. 

**For the case of the Euclidian vacuum the relative signs of the two contributions coming from the expanding 
and contracting PPs can be different depending on whether D is even or odd 9 . But that does not conceptually 
affect our discussion. 
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